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CHAPTER  I 


INTRODUCTION 

Although  there  is  no  general  agreement  on  a  strict  definition  of 
robustness,  most  writers  illicitly  accept  a  qualitative  definition. 

A  general  statement  of  this  definition  is  as  follows:  a  statistical  test 
is  said  to  be  robust  (with  respect  to  an  underlying  assumption  for  some 
class  of  alternatives  to  the  assumption  and  for  a  fixed  sample  size)  if 
the  power  function  of  the  test  (under  any  member  of  the  class  of  alter¬ 
native  ass unpt ions)  is  not  excessively  larger  than  the  power  function  of 
the  test  under  the  original  assumption  for  parameter  values  where  the 
null  hypothesis  is  true  and  not  excessively  smaller  where  the  alternative 
hypothesis  is  true.  A  great  many  tests  are  based  on  an  assumption  of 
normality  of  the  parent  population  and  robustness  with  respect  to  non- 
norrnal  parent  populations  has  received  by  far  the  most  study. 

Classically,  tests  are  formulated  such  that  a  type  I  error,  rejecting 
the  null  hypothesis  when  it  is  true,  is  the  most  critical  error  and  prob¬ 
abilities  of  this  type  of  error  are  strictly  controlled.  Very  often  the 
sample  size  is  controlled  by  physical  considerations  and,  hence,  the 
probability  of  the  other  type  of  error  is  beyond  the  control  of  the 
experimentor.  For  these  reasons,  changes  in  the  probability  of  a  type 
I  error  for  changes  in  the  assumptions  are  of  more  importance  than  other 
points  on  the  power  function  and  most  robustness  studies  have  been  limited 
to  the  null  hypothesis  point  on  the  pcwer  function.  It  should  be  observed 


1 


2 


that  the  central  issue  is  the  distribution  c'  test  statistic  under  the 
alternative  assumptions  being  considered.  The  specific  concern  in  this 
work  is  the  density  of  the  one  sample  t-statistic  without  the  assumption 
of  normality  of  the  parent  population. 

to  annotated  bibliography  of  robustness  studies  in  general  has  teen 
given  by  Grvindarajulu  and  Leslie  [8].  A  survey  of  robustness  studies  of 
the  Student  t-tests,  both  one  sample  and  two  sample,  has  been  given  by 
Hatch  and  Posten  [9] .  We  will  adopt  the  convention  that  references  to 
the  Student  t-statistic  or  test  have  the  underlying  assumption  of  normality 
of  the  parent  population,  while  references  to  the  t-statistic  or  test 
include  no  such  assumption. 

Surprisingly  little  has  been  accomplished  in  deriving  the  exact 
density  or  mass  function  of  the  t-statistic.  Rider  f 15]  derived  the 
density  for  samples  of  size  2  for  a  uniform  population  as  well  as  the 
mass  function  for  various  discrete  uniform  samples  of  size  2,  3,  and  4. 
Perlo  [14]  nas  given  the  density  of  the  t-statistic  for  samples  of  size 
3  from  a  uniform  parent  population.  Geary  [7]  derived  the  t-density  for 
double  exponential  samples  and  Baker  [1]  treated  the  compound  normal  case 
with  equal  variances;  both  for  sample  size  2.  Hotelling  [10]  derived 
the  tails  of  Che  t-density  for  samples  of  size  2  from  a  Cauchy  parent 
peculation.  Ladorman  [11]  derived  the  t-density  for  samples  of  size  2 
from  an  arbitrary  density  wiuh  mean  zero  by  geometric  arguments.  His 
result  is  derived  analytically  here  without  the  assumption  about  the  mean 
but  for  parent  densities  positive  on  the  entire  real  line  and  can  easily 
be  derived  for  the  other  cases  with  results  presented  here. 

Various  approximations  for  the  t-density  have  been  given.  Bartlett 
[2],  Geary  [7],  Gayen  [6]  and  others  have  used  the  first  few  terms  of  an 
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Edgeworth  type  A  series  or  a  Gras-Charlier  series  as  the  density  of  the 
parent  population  and  proceeded  to  derive  the  associated  t-dens-ty.  A 
thorough  account  of  many  such  works  is  given  by  Hatch  and  I*osten  [9j. 
Bradley  [41  worked  in  quite  a  different  way.  He  wrote  the  distribution 
function  of  the  t-statistic  as  an  integral  of  the  joint  density  of  the 
observation  over  the  appropriate  subset  of  Euclidean  n -space,  then  manip¬ 
ulated  the  n-fold  integral.  After  making  sinplifying  assumptions  about 
the  parent  density,  similar  to  those  made  in  this  work,  he  approximated 
the  t-density  with  the  first  few  terms  of  a  series  representation  for 
it.  He  developed  a  computational  technique  and  illustrated  it  with 
Cauchy  and  logistic  parent  populations  for  sample  sixes  of  2,  3,  4  and  5. 

The  general  approach  taken  here  is  to  recursively  derive  the  joint 
density  of  the  sample  mean  and  sum  of  squares  of  deviation  about  the 
sample  mean.  For  n  >  3,  the  recursion  relation  requires  an  integration 
which  is  accomplished  by  application  of  the  mean  value  theorem.  This 
technique  produces  exact  results  for  only  a  certain  class  of  functions 
but  gives  an  approximation  for  others.  This  type  of  application  of  the 
mean  value  theorem  to  carry  out  integrals  promises  to  be  a  powerful 
statistical  technique  with  further  study,  which  is  indicated  in  Appendix  A,. 
A  transformation  from  this  density  to  the  density  of  the  t-statistic  is 
given  here,  along  with  a  symmetry  property  for  the  t-density. 

The  formulas  derived  here  are  illustrated  when  a  member  of  the 
generalized  normal  family  is  the  parent  density.  The  approximation  of 
the  t-density  is  derived  for  these  parent  densities  for  the  case  where 
u  =  0  and  for  all  sample  sizes.  Tables  of  type  I  error  probabilities  are 
given  for  several  specific  members  of  this  family  for  sample  sizes 
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CHAPTER  II 


can  be  verified  directly.  The  superscript  on  T^  and  will  be  suppressed 
where  there  is  no  ambiguity. 

fy.  \  /-.\ 

Suppose  f  (t.,  t.)  is  the  joint  density  of  T,  and  t'  .  X  . , 

n  l  <s  i  i  n+x 

is  independent  of  (X,,  X„,  •••  ,  X  )  and  hence,  is  independent  of 

1  i  n 

and  .  Therefore,  the  joint  density  of  and  Xn+^  is 

f  (t, ,  t.)f(x  .  Using  the  relations  in  (2.1),  we  can  transform 
n  1  2  n+i  * 

[T.^  ,  T^n^ ,  x  ]  -»■  [T^n+^  ,  T^n+^  ,  U]  with  the  auxiliary  variable  being 
l  <s  n+i  l  2 

/n  \ 

defined  by  U  =  T;  -  X  , ,  .  The  inverse  of  this  transformation  is 
1  1  n+1 


¥- 


Ri 


An) 

‘t.  = 

1 

1 

Ai 

T  n+1 

„(n) 

T  .TB 

(n+1) 

n 

2 

2 

n+1 

X  = 

(n+1) 

n 

n+1 

1 

~  n+1 

and  the  jacobian  is  J  *  -1  .  Then  the  density  of  (T*n+^ ,  ,  0]  is 


£n(‘l  +  ST  £2  '  ST  u2)£(tl  -  ST")  , 


<  tx  <  «  ,  t2  >  0  , 


and  the  density  of  [T^n+^ ,  *2n+^  1  is 


fn+l(t 


I'V  *  }s  £n(£l  +  SI  "•  ‘2  '  SI  "2)f(ti  -  SI  u)d"  ' 


<  tj^  <  00  ,  t2  >  0  , 


where  S  is  the  open  interval 

S  =  ^-  ^  (n+1)  t2/n  ,  y(n+l)t2/n 

In  order  to  make  the  range  of  integration  independent  of  t2  and  to 
put  the  recusion  relation  in  a  more  useful  form  we  can  transform 

u  =  ^  <n+l)t2/n  v  . 

. ,  ,  .  ,  . ,  _  _  (n+1)  (n+1)  . 

Then  the  joint  density  of  and  T2  is 
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The  only  changes  required  in  (3.1)  for  parent  densities  that  are  positive 
on  (a,  b)  or  (0,  ®)  are  changes  in  the  limits  on  t;  and  t?  .  For  den¬ 
sities  positive  on  (a,  b) ,  a  <  t^  <  b  and  0  <  /tj  <  J2  mintt^-a,  b-t^)  . 
For  densities  positive  on  (0 ,  ®) ,  0  <  t^  <  ®  and  0  <  /t~  <  *^2  t^  .  These 
results  are  also  proven  by  Craig  [5]. 

Applying  the  relation  (2.2)  we  have 


Let 


then  we  can  rewrite  f^ (t^,  t^)  in  the  more  compact  form 

W  V  *  ^  f  (1-“2>"1/2  ",  £[‘i  *  ^  \3<“>] 

J-l  1=1 


du  .  (3.4) 


If  we  further  restrict  f  to  be  everywhere  continuous,  the  integrand  is 
the  product  of  compositions  of  continuous  functions  on  the  interval 
(-1,  1).  Hence,  the  integrand  is  a  continuous  function  on  (-1,  1)  that 
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is  unbounded  at  the  end  points.  Then  for  each  fixed  (t^,  t 2>  we  can  apply 
the  modified  mean  value  theorem  for  integrals ,  a  statement  and  proof  of 
which  is  given  in  Appendix  A,  to  say  there  exists  (-1,  1)  such  that 


fjitj,  t2>  =  2/3  u  -  ih~1/2  u 


tl  +  ^  ai3<53> 


]• 


(3.5) 


3  2 

where  £  a.  (?)  =  0  and  \  a.,(£.J  -  1  are  easily  verified  identities 


i3  3 


i3  3 


in  gy  as  well  as  the  fact  that  -1  <  a^3  <  1  for  i  =  1,  2,  3  . 

The  existence  of  a  value  ?3  is  guaranteed  for  each  (t^,  t2>  ,  hence, 
C3  =  ^(f,  t^,  t2).  The  part  of  the  integrand  involving  f  in  (3.4)  can 


be  written 


+a“(u’])- 

For  densities  that  are  symmetric  about  zero  and  sufficiently  smooth, 
t-j/Zt”  large  moves  the  range  of  integration  out  in  the  tail  of  the 
density  and  the  multiplicative  factor  tends  to  play  a  smaixer  role  as  the 
ratio  increases.  Then  for  smooth  densities  symmetric  about  zero 
£3  «  C3[f,  |t^//tj  |J  .  Also,  for  densities  that  can  be  written  with 
a  scale  parameter,  the  scale  can  be  made  large  and  the  mass  concentrated 
about  zero  so  that  the  ratio  tj//tJ  is  always  large  and  hence,  £3  «  ^(f) 
since  T  \./Jr2  is  a  scale  invariant  random  variable.  In  any  case,  we  will 
treat  as  a  constant  and  derive  the  values  of  the  coefficients  {a^} 
from  other  considerations.  Then  the  following  results  will  be  exact 
only  for  the  class  of  density  functions  where  £3  is  independent  of  t^ 
and  t2  .  For  the  class  of  densities  where  ?3  is  not  independent  of  t^ 
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and  t 2  results  that  follow  will  yield  approximations  for  the  densities 
of  Cq,  0?2)  and  T  . 

Suppose  we  have  recursively  derived,  for  all  integers  up  to  and 
including  n. 


W  V 


*£*jn  a-q) 


i-4 ' 

2,  2 


n  f(t,  +  ft!  a.  ) 
1  2  m 


-»<t^<<»,t2>o 


(3.6) 


where 


And  for  each  i  <  n  ,  q  is  the  constant  whose  existence  is  guaranteed  by 
the  mean  value  theorem.  Applying  the  recursion  relation  (2.2),  we  have 


-®  <  q  <  ®  ,  t2  >  0  .  (3.8) 
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Again  the  mean  value  theorem  can  be  applied  to  provide  the  existence  of 

C  ^£[-1,  1]  which  is  assumed  to  be  independent  of  t.^  and  t 2,  using  the 

n 

conditions  stated  following  equation  (3.5).  The  identities 

5  2 

and  >  a.  =  1  can  be  verified  directly. 

“  x,n+l 

Then  by  the  strong  principal  of  finite  induction,  the  density  of 

(T^,  T^)  is  given  by  (3.6)  and  the  recursion  relation  for  the  coefficients 

(a.  }  in  terms  of  the  constants  {£.)  is  given  by  (3.7)  for  all  values  of 
xn  x 

n  >  3  .  If  we  take  vacuous  products  to  be  1,  (3.6)  is  valid  for  n  >  2  . 


a.  . .  «  0 
i,n+l 


CHAPTER  IV 


12 


13 


Q  (t)  is  an  even  function  whenever  f  is  an  even  function.  This 
n 

will  be  verified  by  the  fact  that  for  every  n  >  2,  f  (c.,  t_)  -  f  (-t. /  t_) 

n  x  «  n  x  / 

for  all  t^  and  t^,  which  will  first  be  established  by  induction.  Since 
f  is  an  even  function,  the  property  holds  for  n  -  2  from  (3.1).  Supposing 
the  property  holds  for  n,  applying  (2.2) ,  we  have  for  all  t^  and  1 2 

V  =V"n~  fc2  (  fn(fcl  "  uVn(n+l)  '  t2(1_U  ^^l  +  )du  * 


Then  with  the  change  of  variable  u  -u,  we  see  the  property  holds  for 

n+1  and  hence  the  property  holds  for  n  >  2  . 

The  density  of  (T- ,  /FT)  =  (X,  Y)  is  2yf  (x,  y^)  and  we  have,  since 
12  n 


Y  >  0, 


!«  --tx 

I 

0 


2yf  (x,  y  ) dxdy 
n 


Then  using  the  symmetry  property  of  f  and  making  the  change  of  variable 
x  -*■  -x,  we  have 


Pr  y  t  "t  55  Pr  y  >  t  ,  for  all  t  . 


Hence,  T^/ZrJ  is  a  symmetric  random  variable  and  is  an  even  function. 

Comparing  Qn(t)  andQn(-t)  in  (4.1),  one  implication  of  this 
symmetry  might  be  that  {a^n|i  =  1,  2,  •••  ,  n)  =  {-a^n|i  =  1,  2,  *••  ,  n) 
for  some  f's.  When  this  implication  is  not  true,  the  set  {a^}  that  is 
symmetric  about  zero  could  be  used  to  approximate  the  t-density  with 
normalized  to  make  it  a  density  function,  since  the  symmetric  set  of 
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coefficients  does  make  an  even  function.  This  approximation  could 
also  be  applied  to  parent  densi.  as  that  are  not  even  functions.  The 
coefficient  set  {a^}  that  is  synnvetric  is  dispersed  on  the  interval 
(-1,  1)  and  hence  would  be  a  reasonable  approximation  of  the  true  coeffi¬ 
cient  set  in  the  absence  of  additional  information  about  the  true  set. 
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DETERMINATION  OF  SYMMETRIC  [a^j 


By  the  relation  (3.7) ,  we  can  see  that  the  coefficient  set 
{aji  =  1,  2,  •••  ,  n}  has  the  property  a^-a^  =  (a.  ^ 
for  i,  j  =  1,  2,  ,  n-1  and  hence  any  ordering  of  {ai  n-1| i  =  1,  2, 

•••  ,  n-1}  must  give  the  same  ordering  of  { ain I  i  =  1#  2,  n~l>  •  When 
we  have  solved  for  £  .  and  {a.  .}  we  will  order  {a.  .}  and  then  we 

will  have  only  three  considerations  to  determine  the  position  of  ann  in 
the  order.  For  n  =  2k+l,  the  (k+l)st  value  in  the  completely  ordered 
set  must  be  zero,  i.e.,  a^  n  =  0,  n  =  0  or  a^n  =  0  •  For  n  =  2k, 

the  extreme  values  in  the  completely  ordered  set  must  differ  exactly  in 


sign,  i.e.,  a^  =  -a^^  ,  =  -ann  ,  or  a 


—  **s  .  • 

nn  n-l,n 


For  n  =  3,  the  coefficients  are  given,  as  functions  of  by  (3.3) 

and  the  conditions  above  yield  the  solutions  =  0,  -/3/2,  /J/2  .  All 

three  values  of  yield  the  coefficient  set  (a  }  =  {l//2,  0,  -1//2  ]  . 

3  J*  *3 

2 

We  can  discriminate  between  =  ^  or  3/4  by  requiring  Q^ft)  h°  ke  a 

density  function.  For  n  =  4,  the  coefficients  are  given  by  (3.7)  and  the 

conditions  above  yield  5^  =  0,  -*^3/5  ,  /3/5  .  The  values  =  ^3/5  ,  -/3/5 

yield  the  coefficient  set  j  3//20*  ,  1/^20  ,  -I//20"  ,  -3/^20  j  and  the 

value  =  0  yields  the  coefficient  set  {  1//2  ,  0,  0,  -I//2J  .  We  will 

assume  that  %  /  0,  n  >  4,  and  handle  the  other  cases  later.  By  inspec- 
n  ~ 

tion  of  these  cases,  we  can  set  up  the  induction  hypotheses 
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THE  GENERALIZED  KORMAL  DISTRIBUTION 

The  generalized  normal  distribution  has  been  considered  by  various 
writers  in  connection  with  robustness  studies.  The  density  is 


3+8 

,  2  _ 

(3+6) 

-1 

1|  x-u  i1^] 

2  r 

It r 

exp 

r  2  \—\ 

[ 

«*.•.«  r|  2Tj  'M*’-f|TL| 

’*•<*<»,  -l<Q<lt-?»<v<~rO>0 


Due  to  the  scale  invariance  of  T,  we  will  consider  a  -  1,  without  loss  of 
generality.  We  will  further  consider  the  null  case  where  y  =  0,  and  f  is 
an  even  function.  Applying  (4-1)  we  have 

\  n  |  „  I  mHP?) 


where 


n 

t 

2 

1+B 

i 

1 

/n(n-l)  in 

k 

- 

d+e)r[nJ^-)J  n  (i-c?) 

c  (6)  = - - - 

rn(^) 

Since  must  be  an  even  function,  we  also  have 


(6.1) 


\  n 
/  V 

t 

2 

1+8 

-m 

1  i 
/  1 

A (n-1)  in 

17 


18 


19 


n  3 

of  all  odd  powers  of  t  cost  be  zero.  The  solutions  of  J  a .7  *  0  are 

,  iy> 

i 


Z  =0,  +  w— 2-  which  are  fee  values  that  produce  the  syapetric  constants 

IS  •  I  Sri 

discussed  in  the  previous  chapter.  Clearly,  these  syefcric  constants  have 
5  2k+l 

the  property  that  )  a.  =0  for  k  *  0,  1,  2,  •••  .  t&en  2/ (1+8)  *  2k+l, 

1  in 

we  can  cccpare  (6.1)  and  (6.2}  to  see  that 


n 

£ 


tr.  (n-i) 


+  a. 
in 


2k*l  n 


/a  (n*l) 


-  a. 


in 


2k+3. 


(6.3) 


is  an  identity  in  t  .  For  k  =  0,  we  can  90  through  a  tedious  process  of 

letting  t//n(n-l)  take  values  between  tbs  a.  *s  to  show  that  the  set  of 

in 

coefficients  {a.  }  cast  be  symmetric.  For  k  >  1,  we  can  let 
in  - 

t/»n  (n-1)  >  max  j  |  a.  j  J  and  expand  the  tens  by  the  binomial  tneoree-  Ccm- 


1* 


paring  terns  in  (6.3)  we  can  again  see  that  \  a^*T*  *  0  for  k  =  1,  2,***, 

x 

and  hence  that  the  symmetric  set  eaist  be  the  proper  cne.  Than  the 

symmetric  set  of  constants  must  apply  for  all  3  /  0  such  that  2/ (1+8)  is 


an  integer. 

The  set  of  symmetric  coefficients  was  used  in  (4.1)  for  the  values 
of  8  where  2/ (1+8)  =  1,  4/3,  3/2,  7/3,  5/2,  4,  16  .  The  appropriate  values 
of  =  0,  3/4  where  used  to  numerically  integrate  Q^(t|S) ,  since  this 
must  be  a  density  function,  neither  of  the  values  gave  a  density  function. 
This  shows  that  f{xjo  ’  3)  is  not  in  the  class  of  parent  densities  for 


which  (4.1)  is  the  exact  t-density,  but  (4.1)  does  represent  an  approxi¬ 


mation  for  this  density. 
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The  set  {£  )  ,  with  the  associated  set  of  coefficients  (a.  }  ,  was 
n  m 

fitted  sequentially  for  the  values  of  3  given  above  in  such  a  way  as  to 

sake  the  integral  of  Q^CtjS)  closest  to  i  .  The  sanple  sizes  considered 

were  n  *  2,  3,  •••  ,  31  and/  of  course,  for  n  =  2  the  density  of  Q2(tj&) 

2 

is  exact  for  all  3  .  The  set  {£  }  used  was  £  -  3/(n+l)  except  for 

n  n 

the  0's  where  2/(l+Sil  <  2  and  in  this  case,  =  S13  =  Cig  =  Cg  =  0 

were  the  only  changes. 

The  tables  that  follow  are  tables  of 


e(n,  3) 


(6.4) 


where  t  (n-1)  is  the  critical  point  of  the  Student  t-density  for  a  normal 
a 

parent  population.  The  ta\n~l)  were  taken  from  [12],  for  the  five  decimal 
place  accuracy,  where  available  and  from  [13],  otherwise,  a  is  interpreted 
as  the  true  probability  of  a  type  I  error  under  the  approximation  Qn  (t  |  B ) 
when  a  was  the  advertised  probability  and  the  false  assumption  of  normality 
was  utilized.  The  compute r  program  that  was  used  to  carry  out  the  cal¬ 
culations  is  listed  and  explained  in  Appendix  B. 
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TABLE  1 


n-1 

a=.100 

a=.050 

a=.025 

a=.010 

ct=.005 

1 

,08129 

.03962 

.01969 

.00736 

.00393 

2 

.08885 

.03705 

.01707 

.00651 

.00321 

3 

,10313 

.04453 

.01801 

.00620 

.00291 

4 

.11142 

.05186 

.02237 

,00675 

.00296 

5 

.11945 

.05812 

.02620 

.00831 

.00336 

6 

.09308 

.04741 

.02211 

.00778 

.00341 

7 

.11259 

.05742 

.02952 

.01143 

.00497 

8 

. 12572 

.07042 

.03744 

.01439 

.00647 

9 

.09638 

.05148 

.02760 

.01103 

.00529 

10 

.11456 

.06460 

.03424 

.01440 

.00698 

11 

.13114 

.07277 

.04011 

.01692 

.00872 

12 

.10234 

.05411 

.02889 

.01214 

.00615 

13 

.11645 

.06465 

.03498 

.01544 

.00827 

14 

. 12540 

.06987 

.03948 

.01328 

.00966 

15 

.13046 

.07618 

.04405 

.02012 

.01085 

16 

.14016 

.08252 

.04719 

.02175 

.01185 

17 

.  14570 

.08531 

.04888 

.02301 

.01268 

18 

.11393 

.06603 

.03672 

.01651 

.00904 

19 

.12571 

.07091 

.03978 

.01857 

.01003 

20 

.12967 

.07412 

.04259 

.01961 

.01069 

21 

,13  223 

.07695 

.04360 

.02015 

.01112 

22 

.  13677 

.07853 

.04467 

.02102 

.01169 

23 

.13763 

.07917 

.04562 

.02161 

.01219 

24 

.13855 

.03067 

.04667 

.02250 

.01273 

25 

.13911 

.08134 

.04762 

.02304 

.01322 

26 

. 14047 

.08307 

.04870 

.02395 

.01372 

27 

. 14290 

.08472 

.05021 

.02466 

.01429 

28 

.11954 

.06776 

,03850 

.01822 

.01027 

29 

.12415 

.07087 

.04079 

.01947 

.01105 

30 

.12768 

.07392 

.04263 

.02048 

.01155 

2 


TABLE  2 

2  4 

1+B  3 


a=.100 

a=.G50 

a=.025 

a=.010 

ct=.005  1 
1 

.08936 

.04395 

.02188 

.00874 

.00437 i 

.09323 

.04220 

.02006 

.00778 

.00385  1 

.10125 

.04637 

.02074 

.00759 

.00365 

.10652 

.05066 

.02314 

.00799 

.00370 

. 11082 

.05425 

.02542 

.00390 

.00398 

.09771 

.04366 

.02330 

.00855 

.00293 

.10771 

.05458 

.02734 

.01049 

.00483 

.11485 

.06063 

.03100 

.01199 

.00562 

. 10044 

.05165 

.02646 

.01045 

.00504 

. 10926 

.05801 

.02994 

.01214 

.00593 

.11717 

.06237 

.03276 

.01344 

.00673 

.10336 

.05355 

.02766 

.01130 

.00565 

.11010 

.05339 

.03053 

.01233 

.00657 

.11496 

.06153 

.03288 

.01412 

.00725 

.11833 

.06462 

.03496 

,01505 

.00730 

.12260 

.06741 

.03647 

.01580 

.00825 

. 12543 

.06399 

.03743 

.01637 

.00860 

.11317 

.06049 

.03219 

.01373 

.00719 

.11601 

.06251 

.03351 

.01457 

.00761 

.11306 

.06409 

.03471 

.015 10 

.00794 

.11965 

.06548 

.03543 

.01549 

.00319 

.12157 

.06644 

.03605 

.01590 

.00846 

.12248 

.06709 

.0366  1 

.01624 

.00870 

. 12322 

.06783 

.03716 

.01661 

.00393 

.12332 

.06339 

.03766 

.01690 

.00914 

. 12460 

.06914 

.03817 

.01724 

.00934 

.  12565 

.06990 

.03S74 

.01754 

.00954 

.11541 

.06259 

.03393 

.01502 

.00305 

.11691 

.06368 

.03472 

.01545 

.00332 

.11330 

.06431 

.03543 

.01583 

.00352 

23 


TABLE  3 

2  _  2 
1+8  "  2 


n-1 

- -  — . 

ct^.100 

a=.050 

a=.025 

a=.010 

a=.005 

1 

.09258 

.04573 

.02279 

.00911 

.00455 

2 

.09512 

.04442 

.02141 

.00838 

.00416 

3 

. 10070 

.04732 

.02194 

.00824 

.00400 

4 

. 10450 

.05033 

.02360 

.00855 

.00405 

5 

. 10748 

.05284 

.02520 

.00919 

.00426 

6 

.09896 

.04913 

.02379 

.00893 

.00420 

7 

. 10559 

.05327 

.02655 

.01025 

.00484 

8 

. 11044 

.05716 

.02890 

.01125 

.00537 

o 

✓ 

. 10 107 

.05142 

.02602 

.01027 

.00499 

10 

.10633 

.05556 

.02835 

.01139 

.00559 

11 

.11191 

.05847 

.03019 

.01226 

.00611 

12 

.10304 

.05284 

.02699 

.01092 

.00544 

13 

.10741 

.05595 

.02889 

.01191 

.00603 

14 

.11070 

.05813 

.03041 

.01273 

.00647 

15 

.11311 

.06014 

.03173 

.01333 

.00632 

16 

.11573 

.06139 

.03270 

.01382 

.00710 

17 

.11763 

.06295 

.03335 

.01413 

.00732 

13 

. 10994 

.05770 

.03016 

.01264 

.00650 

19 

.11169 

.05396 

.03093 

.01311 

.00676 

20 

.11305 

.05999 

.03172 

.01347 

.00697 

21 

.11416 

.06039 

.03223 

.01374 

.00714 

22 

.11536 

.06156 

.03267 

.01400 

.00731 

23 

.11605 

.06205 

.03305 

.01423 

.00746 

24 

.11661 

.06255 

.03341 

.01446 

.00761 

25 

. 11709 

.06297 

.03375 

.01465 

.00774 

26 

.11763 

.06344 

.03407 

.01486 

.00736 

27 

.11328 

.06391 

.03441 

.01503 

.00797 

23 

.11200 

.05943 

.03155 

.01357 

.00713 

29 

.11231 

.06009 

.03199 

.01331 

.00728 

30 

. 11362 

.06074 

.03240 

.0 1404 

.00740 

24 


TABLE  4 

2  2 
1+g  “  3 


n-1 

.  1  ■'  ■  ' 

a=.100 

a= . 050 

ct=.025 

a=.010 

a=.005 

1 

.  10361 

.05220 

.026 15 

.01047 

. 00523 

2 

. 10268 

.05298 

.02703 

.01096 

.00551 

3 

.10003 

.05160 

.02668 

.01103 

.00561 

4 

.09802 

.05012 

.02588 

.01082 

.00556 

5 

.09649 

.04888 

.02508 

.0104  3 

.00542 

6 

.09533 

.04791 

.02440 

.01014 

.00524 

7 

.09450 

.04715 

.02384 

.00983 

.00507 

8 

.09383 

.04653 

.02338 

.00957 

.00491 

9 

.09326 

.04603 

.02200 

.00934 

.00477 

10 

.09231 

.04562 

.02269 

.00915 

.00465 

11 

.09244 

.04528 

.02242 

.00399 

.00455 

12 

.09212 

.04498 

.02219 

.00385 

.00446 

13 

.09184 

.04473 

.02199 

.00873 

.00433 

14 

.0916  3 

.04451 

.02132 

.00363 

.00431 

15 

.09143 

.04432 

.02167 

.00353 

.00425 

16 

.09121 

.04415 

.02154 

.00845 

.00420 

17 

.09106 

.04400 

.02142 

.00333 

.00415 

18 

.09091 

.04388 

.02132 

.00331 

.00411 

19 

.09079 

.04375 

.02122 

.00325 

.00407 

20 

.09068 

.04364 

.02114 

.00320 

.00403 

21 

.09056 

.04354 

.02106 

.00815 

.00400 

22 

.09047 

.04345 

.02099 

.00310 

.00397 

23 

.09037 

.04336 

.02092 

.00806 

.00394 

24 

.09030 

.04329 

.02036 

.00302 

.00392 

25 

.09023 

.04322 

.02031 

.00799 

.00390 

26 

.09015 

.04315 

.02075 

.00795 

.00333 

27 

.09009 

.04309 

.02071 

.00792 

.00386 

28 

.09003 

.04304 

.02066 

, 007  89 

.00333 

29 

.08993 

.04293 

.02062 

.00737 

.00382 

30 

.03992 

.04294 

.02053 

.00734 

.003S0 

25 


TABLE  5 

2  5 

1+6  =  2 


n-1 

o=.100 

a=.050 

ct=.025 

a=.010 

a=.005  | 

1 

.10512 

.05315 

.02666 

.01067 

.00534 

2 

.103  36 

.05428 

.02795 

.01140 

.00574 

3 _ 

.10013 

.05233 

.02743 

.01150 

.00539 

4 

.09724 

.05023 

.02629 

.01119 

.00582 

5 

.09505 

.04847 

.02515 

.01070 

.00561 

6 

.09343 

.04708 

.02418 

.01021 

.00536 

7 

.09215 

.04599 

.02339 

.00973 

.00511 

3 

.09117 

.04510 

.02274 

.00941 

.00438 

9 

.09035 

.04439 

.02220 

.00909 

.00469 

10 

.03970 

.04379 

.02175 

.00382 

.00452 

11 

.089 14 

.04330 

.02137 

.00360 

.00433 

12 

.08368 

.04287 

.02105 

.00340 

.00425 

13 

.08827 

.04251 

.02077 

.00823 

.00415 

14 

.08793 

.04219 

.02053 

.00308 

.00405 

15 

.08762 

.04192 

.02032 

.00795 

.00397 

16 

.08735 

.04  lt>  8 

.02013 

.00733 

.00339 

17 

.08712 

.04146 

.01996 

.00774 

.00333 

IS 

.03690 

.04127 

.01931 

.00764 

.00377 

19 

.08672 

.04109 

.01963 

.00756 

.00371 

20 

.03655 

.04093 

.01956 

.00749 

.00366 

21 

.03638 

.04079 

.01945 

.00742 

.00362 

22 

.03625 

.04066 

.01934 

.00735 

.00358 

23 

.0  3611 

.04054 

.01925 

.00730 

.00354 

24 

.03600 

.04043 

.01917 

.00724 

.00351 

25 

.08589 

.04033 

.01909 

.00720 

.00348 

26 

.08577 

.04023 

.01901 

.00715 

.00345 

27 

.08563 

.04014 

.01895 

.00711 

.00342 

28 

,08560 

.04007 

.01888 

.00707 

.00340 

29 

.08552 

.03999 

.01382 

.00703 

.00337 

30 

.08544 

.03^92 

.01877 

.00700 

.00335 

26 


TABLE  6 


n-1 

a=. 100 

a=.050 

a=.025 

a=.010 

a=.CG5 

1 

.11349 

.05398 

.02983 

.01197 

.00599 

2 

.11074 

.06206 

.03332 

.01447 

.00743 

3 

.10170 

.05718. 

.03222 

.01472 

.00794 

4 

.09334 

.05170 

.029 12 

.01365 

.G0760 

5 

.08763 

.04701 

.02612 

.01225 

.00691 

6 

.03275 

.04321 

.02356 

.01093 

.00617 

7 

.07888 

.04015 

.02144 

.00977 

.00549 

8 

.07575 

.03765 

.01970 

.00880 

.00489 

9 

.07318 

.03560 

.01827 

.00799 

.00439 

10 

.07104 

.03390 

.01707 

.00731 

.00397 

11 

.06925 

.03246 

.01606 

.00674 

.00361 

12 

.06771 

.03124 

.01521 

.00625 

.00331 

13 

.06633 

.03018 

.01447 

.00534 

.00305 

14 

,06523 

.02927 

.01334 

.00548 

.00283 

15 

.06421 

.02347 

.01323 

.00518 

.00264 

16 

.06331 

.02776 

.01280 

.00491 

.00247 

17 

.06252 

.02713 

.01237 

.00467 

.00233 

18 

.06180 

.02657 

.01198 

.00446 

.00220 

19 

.06116 

.02607 

.01164 

.00427 

.00208 

20 

.06053 

.02561 

.01133 

.00411 

.00198 

21 

.06005 

.02520 

.01105 

.00396 

.00189 

22 

.05957 

.02432 

.01080 

.00382 

,00131 

23 

.05912 

.02448 

.01057 

.00370 

.00174 

24 

.05372 

.024 16 

.01036 

.00358 

.00167 

25 

.05834 

.02387 

.01016 

.00343 

.00161 

26 

.05799 

.02360 

.00993 

.00339 

.00156 

27 

.05767 

.02335 

.00982 

.00330 

.00150 

28 

.05737 

.02312 

.00966 

.00322 

.00146 

29 

.05709 

.02290 

.00952 

.00315 

.00142 

30 

.05632 

.02270 

.00939 

.00308 

.00133 

27 


TABLE  7 


n-1 

a= . 100 

a=.050 

a=.025 

a=.0!0 

a=.005 

1 

.12187 

.06730 

.03534 

.01442 

.00724 

2 

.11510 

.06936 

.04097 

.01946 

.01065 

3 

.  10222 

.06194 

.03781 

.01935 

.01139 

4 

.09037 

.05386 

.03290 

.01724 

.01048 

5 

.03024 

.04663 

.02819 

.01479 

.00909 

6 

.07167 

.04059 

.02410 

.01252 

.00771 

7 

.06440 

.03547 

.02066 

.01057 

.00647 

8 

.05820 

.03117 

.01778 

.00893 

.00542 

9 

.052S7 

.02754 

.01538 

.00757 

.00455 

10 

.04  ^27 

.02446 

.01338 

.00645 

.00383 

11 

.04428 

.02184 

.01169 

.00552 

.00324 

12 

.04079 

.01960 

.01027 

.00474 

.00275 

13 

.03771 

.01766 

.00907 

.00410 

.00234 

14 

.03501 

.01599 

.00804 

.00355 

.00201 

15 

.03260 

.01453 

.00716 

.003 10 

.00172 

16 

.03045 

.01325 

.00640 

.00271 

.00149 

17 

.02353 

.01212 

.00575 

.00238 

.00129 

18 

.02631 

.01113 

.00518 

.00210 

.00112 

19 

.02525 

.01025 

.00468 

.00186 

.00098 

20 

.02384 

.00947 

.00424 

.00165 

.00086 

21 

.02256 

.00876 

.00386 

.00147 

.00076 

22 

.02139 

.00814 

.00352 

.00132 

.00067 

23 

.02032 

.00757 

.00322 

.00118 

.00059 

24 

.01935 

.00706 

.00295 

.00106 

.00053 

25 

.01845 

.00660 

.00271 

.00096 

.00047 

26 

.01762 

.006  IS 

.00249 

.00086 

.00042 

27 

.01635 

.00580 

.00230 

.00073 

.00037 

23 

.01615 

.00545 

.00213 

.00071 

.00034 

29 

.01549 

.00513 

.00197 

.00065 

.00030 

30 

.01438 

.00484 

.00183 

.00059 

.00027 

-CHAPTER  VII 


CCWCLDSIGKS 

As  was  previously  ntu.  -cnc-dt  many  authcrs  have  concluded  that  kur- 
tosis  is  of  pt -unary  Importance  in  the  robustness  of  the  Student  t-tesi. 

With  our  appr rwrinati  on ,  we  mi  ht  extend  this  to  say  that  for  parent  pop¬ 
ulations  w'th  kurtosis  >  3,  which  is  the  normal  density  value,  the  tests 
axe  optimisric  in  the  type  I  err, or  probabilities;  that  is,  the  real  a 
value  is  larger  than  the  advertised  value.  For  parent  populations  with 
kurlosis  <  3  the  tests  are  conservative ;  that  is,  the  real  a  value  is 
smaller  than  the  advertised  value. 

Box  and  Tiao  [3]  concluded  that  the  Student  t-test  is  a  good  approxi¬ 
mation  for  all  members  of  the  generalized  normal  family.  They  used  Bayesian 
techniques  with  specific  families  of  prior  densities  for  y  and  &  in  their 
work.  The  present  work  seems  to  reinforce  their  conclusions  since  the  same 
conclusions  follow  from  very  different  approaches. 

It  should  be  observed  that  the  conclusions  stated  here  are  based  on 
an  approximation  of  the  t-density  for  a  particular  family  of  parent  den¬ 
sities.  Further,  the  approximation  was  thoroughly  investigated  for  only 
a  few  members  of  this  family  and  nc  estimate  is  given  of  the  precision 
of  the  approximation. 


20 


A2?£30IX  n 


Far  th?  sake  of  completeness,  the 


tiseoua  in  its  stndud 


for*  is t  stated.  'She  proof  is  onitted  since  it  can  be 
elementary  calculus  textbooks. 


Theorem;  If  g  is  a  ccctinocat  fuctioc  on  the  interval  ga,  b]>  then 

ib 

there  exists  Ceja,  b!  such  that  1  g(x)dx  *  Cb-a)g{£)  - 

^a 

She  slightly  Modified  version  of  this  theorem  that  is  applied  for  equation 
(3.5)  is  as  follows. 


Theory:  If  g  is  a  continuous  function  on  the  interval  (a,  b)  , 

i  *  i  -  f ■ b 

g(x)  =  .  g(x)  *  •  and  1  g(x)dx  <  *  ,  then  there  exists 

x  -*  a  x-»b  ’  . 

fb  'a 

C,»  C2e  b)  such  that  j  g(x)dx  *  (b-aJgiC^J  *  0>-a)g(C2)  i*ere 

*  r 

H  *  v2  * 

Proof:  Since  j  is  continuous  on  (a,  b)  and  positively  unbounded  at  the 
endpoints ,  g  is  bounded  below  on  {a,  b).  Bence,  g  has  an  absolute 
saninun  cn  {a,  b) ,  that  is ,  there  exists  m  (a,  b)  such  that 
g(x)  >  g(n)  for  all  xe  (a,  b).  Then 


g(x)dx 


(A.l) 


Suppose  equality  holds  in  (A. 1)  and  consider  the  function 
h  (x)  =  g(x)  -  i,(n/ ,  which  is  bounded  below  by  zero  and  continuous 
on  (a,  b)  and  positively  unbounded  at  the  end  points.  We  can 
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aggESEHI  B 


The  ccngisggr  program  used  to  calculate  the  probabilities  given  I® 

the  tables  vas  roa  oa  a  SRP3C  11C8-  ffo r  each  @  rrwriiirrrfT,  the  roping 

tine  was  slisfitly  under  1  aaaate.  She  program  given  bare  is  sot  cf&Saml 

for  the  cachic-e ,  bet  it  is  is  a  form  general  enougfi  to  accomodate  asst 

consideraticas  of  the  generalized  nnzaai  drstribstiao  with  p  *  0  - 

The  traosforaatice  x  *  k/{1+o)  was  made  to  sake  toe  range  of 

integration  finite  where  o  is  the  original  variable  of  integration.  'She 

function  given  in  ,,6.1)  was  factored  slightly  and  the  coefficient* 

b.  =  a.  mfe-ij  were  used, 
in  m 


DOUBLE  PEECISIOK  T,  BSSA,  BETAS,  P2B,  GAKKZX,  3,  E,  CP.  SB, 
IGLNi,  G£K2,  C,  PCS',  AKS,  JL 
DI3EKSLC?*  TOSO)  SETA(  )  PSOS(  1  E(31) 

OOJSS3S  SETA,  S,  3(31) 

EXTEUfiAL  EC? 

DKA  T/(ta(n-l)  values) 

DATA  BETAS/ (S  values  to  be  considered) 

DATA  2/(5n  values,  n=2,***,31  with 
DO  5»=1,  (no.  of  S  values) 

XIP=1.QD0 

K=0 

BSTA^BEATAS  (H) 

B(J)=1.CD0 

B(2>=-B(l> 

DO  421=2,  31 
DK=CBI£{?*) 

IF{K.LE.2>  GO  TO  2 

J=N-i 

DO  1  J~1,J 

B  (I)  =E  (S)  +B  (I)  ‘DSQKT  C  (DM*  (1-E  (N)  *E  fti)  )  )  /(DH-2)  ) 

1  CONTIKOE 

B(N)=-(Dft-l)*(2  (H} 


31 


32 


2  GSS1«GMC£ES  (DJI*  (1+323*) /2J 
C!A2*-£K*3XX?jl^  (  (2+323A)  /2J 

33?*XI?*£EXP({  0#-2)/4)  «a£CG(l-2£n>  *2G®)  )  ) 

C*XSF*C  (1+223-A)  /<)  *5£3?  ' (2K/2)  *ELOS  (51)  +  {  (0B-1J/2) 
1*2MG  cES-l;  -K2Jll-KZJ£2 
CXL2,  gCS©  48CFCT,  i.OSO,  0.030,  MIS) 
?3E3{!)«2*C»WS 

HX>  3  >?,  (so.  of  a's  being  considered) 

BW  CEfrj-ll/fff  tt+J-I) +li 

012,  {*£M>  4aCFCT,  l.CDC,  03,  MSS) 

F1!0S£3J»C*35SS/?a3S{l) 

3  cossms 

waraece,  iosojs,  cp°gb(J),  j«i,  fan.  of  o's  +D) 
2000  SOBXKTCSx,  X4,  10*,  ‘no.  off  a*s  +1)210.5,/) 
E=£*teo.  off  o's) 

4  C0BEE55S 

5  Q3S315EE 
StQP 
230 


FEBOTICE  fCT{0) 

000512  PKH^SICEf  3,  X,  *,  FC2,  BETA,  B 
G08SC3K  BETA,  3,  B(32) 

rar=o.Oi» 

EO  1  1=1,* 

WJ*as(x*B(i)) 

5C2=2£3?(  (2/(  1+BEEA)  )  *VWG(t)  )+FCT 
1  C0EFIE5E 

?C?=ESrP(-{3S*  (l+aETA)/2)  *KX3G(2CT)  )/C  (1-0)  *(1-3)  ) 

PEZCZS 

ESS 


JEfflO  48{?CT,  CL,  LL,  ASS)  is  s  48  point  Gauss- Legendre  integration 
subroutine  where  FCX  carries  tile  function  values,  tJL  is  the  upper  Unit 
of  integration,  LL  is  the  lower  limit  of  integration  and  MIS  returns  the 
value  of  the  integral.  GaKAlS(X)  is  a  function  that  calculates  the  natura 
logarithm  of  the  garasa  function  of  x,  using  Bernoulli  nuabers. 
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